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The  two-dimensional,  Eulerian  computer  code  HELP 
I (^drodynamics  ELastic  Plastic)  developed  at  Systems,  Science 

i and  Software  (S*)  has  been  modified  to  allow  for  the  inclusion 

iof  an  explicit  reaction  rate  for  modeling  partially  reacted 
states  of  high  explosives.  Calculations  of  the  detonation 
I thresholds  and  velocity  losses  of  steel  spheres  penetrating 

] slabs  of  Composition  B and  PBX-9404  have  been  carried  out. 

Empirical  reaction-rate  models  based  on  experimental,  one- 
dimensional run-to-detonation  data  (Reference  1)  give  good 
i correlations  with  the  experimental  thresholds  and  velocity- 

I loss  data  as  obtained  at  the  Ballistic  Research  Laboratories 

j (BRL)  . 

The  new  capability  provided  by  the  code  should  be  of 
I use  in  several  problem  areas  of  interest  to  the  Army.  In 

' particular,  work  is  already  underway  to  determine  the  effect 

of  explosive  reaction  rates  on  shaped-charge  jet  formation. 


Section  II  of  this  report  contains  a description  of 
the  reaction-rate  models  and  the  correlations  of  calculated 
and  experimental  velocity-loss  and  detonation  threshold 
data.  The  equation  of  state  and  constitutive  models  used  in 
the  calculations  are  described  in  Section  III.  The  reaction- 
rate  models  and  their  experimental  basis  are  discussed  in 
Section  IV. 
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II.  HELP  CODE  CALCULATIONS 

The  two-dimensional,  Eulerian  HELP  code  (Reference  2) 
has  been  employed  extensively  in  a variety  of  studies  to 
calculate  the  material  flow  associated  with  both  the  impact 
process  and  the  detonation  of  high  explosives  as  well  as  the 
interaction  between  detonated  products  and  inert  material. 
However,  prior  to  the  present  work,  the  explosive  computational 
cell  was  treated  as  either  completely  unreacted  or  completely 
reacted.  In  addition,  it  had  to  be  decided  before  hand 
whether  a given  cell  was  to  detonate,  and  if  so,  at  what 
time.  The  modification  to  allow  for  the  use  of  an  explicit 
reaction  rate  enables  one  to  treat  partially  reacted  states 
of  the  HE.  Correlations  can  be  made  of  the  degree  of  decomposi- 
tion with  the  impact  velocity  of  a projectile,  the  initial  pro- 
jectile/target configuration  and  the  material  properties  of 
both  the  HE  and  the  projectile.  In  particular,  the  use  of  a 
realistic  rate  law  enables  one  to  determine  whether  a given 
impact  and  configuration  will  lead  to  only  partial  decomposi- 
tion or  to  complete  detonation. 

2.1  CORRELATIONS  WITH  EXPERIMENT 

In  the  present  work  the  modified  HELP  code  was  used 
for  a series  of  calculations  consisting  of  a 1.3  cm  diameter 
steel  sphere  striking  at  normal  incidence  a 10  cm  diameter, 

2.5  cm  thick  cylindrical  slab  of  Composition  B or  PBX-9404. 

The  initial  calculational  configuration  is  shown  in  Figure  1. 

The  quantities  of  particular  interest  are  the  velocity  loss 
of  the  projectile,  os  a function  of  the  initial 

striking  velocity,  V , and  the  minimum  striking  velocity  re- 
quired  to  produce  complete  detonation  of  the  slab.  A compari- 
son of  calculated  velocity  losses  with  experimental  data  as 
supplied  by  the  BRL  (Reference  3)  is  given  in  Figure  2. 

Several  of  the  calculations,  as  indicated  in  the  figure,  were 
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HIGH  EXPLOSIVE 


Figure  1.  Initial  configuration  for  the  numerical  calcula 
tions  involving  inert  spheres  perforating  high 
explosive  targets. 


Reactive  EOS,  Detonated 


Velocity  loss  versus  striking  velocity  for  a 1.3  cm  diameter  steel  sphere 
impacting  a 10  cm  diameter,  2.5  cm  thick  cylinder  of  high  explosive. 


performed  using  the  completely  unreacted  equation  of  state 
for  the  HE,  as  given  in  Section  3.2.  The  agreement  of  these 
points  with  the  experimental  curves  is  completely  satisfactory. 
Four  of  the  points  in  Figure  2 are  the  result  of  calculations 
using  the  reactive  equation  of  state  for  HE  (Section  3.2)  and 
the  reaction  rate  models  as  given  in  Figure  3,  which  is  a plot 
of  the  fractional  rate  of  decomposition  as  a function  of  pres- 
sure. As  can  be  seen  from  Figure  2,  the  results  of  these  cal- 
culations are  also  in  good  agreement  with  the  experimental 
results. 

A summary  of  the  results  of  the  calculations  for  PBX- 
9404  is  given  in  Table  I.  The  reaction  rates  used  in  these 
calculations  are  those  as  given  by  the  "Forest  Fire"  model, 
proposed  by  Mader  and  Forest  (Reference  1) . A brief  descrip- 
tion of  the  model  is  given  in  Section  4.1.  The  reaction 
rates  are  obtained  from  experimental,  one-dimensional  run-to- 
detonation  data  and  are  a function  of  a single  equation  of 
state  variable  which  may  be  taken  either  as  the  pressure  or 
the  temperature.  In  the  present  calculation  pressure  was 
used  as  the  independent  variable.  The  Forest  Fire  reaction 
rates  for  PBX-9404  are  plotted  as  a function  of  pressure  in 
Figure  3.  The  use  of  this  model  in  the  HELP  code  calculations 
results  in  a transition  from  partial  to  complete  reaction  for 
striking  velocities  in  the  range  1.0  to  1.1  km/sec,  in  good 
agreement  with  the  experimental  results. 

A summary  of  the  results  of  calculations  for  Composition 
B is  given  in  Table  II.  For  Composition  B,  the  use  of  the 
Forest  Fire  reaction  rates  results  in  only  partial  decomposi- 
tion at  a striking  velocity  of  2.0  km/sec.  Increasing  the 
reaction  rates  by  a factor  of  1.5  at  each  pressure  brings 
the  calculated  detonation  threshold  to  between  1.7  and  1.8 
km/sec,  in  agreement  with  the  experimental  results.  These 
modified  reaction  rates  for  Composition  B are  plotted  in 
Figure  3. 
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TABLE  I 

HELP  CODE  RESULTS  FOR  THE  PENETRATION  OF  PBX-9404 
BY  A 1.3  CM  DIAMETER  STEEL  SPHERE 


R-3512 


TABLE  II 

HELP  CODE  RESULTS  FOR  THE  PENETRATION  OF  COMPOSITION  B 
BY  A 1.3  CM  DIAMETER  STEEL  SPHERE 


REACTION 

RATE  MODEL 

Inert 

Forest  Fire 

FF  X 

1.5 

FF  X 

1.5 

FF  X 

1.5 

FF  X 

1.5 

FF  X 

1.5 

FF  X 

1.5 

X 10 

FF  X 

1.5 

Y 

= 0 

s 

INITIAL 

VELOCITY 

(km/sec) 


VELOCITY 

LOSS 

(km/sec) 


0.56 


.0 

0.85 

0.39 

0.28* 

0.28* 


DETONATION? 


Velocity  loss  at  10  ysec 
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In  summary,  the  equations  of  state  and  constitutive 
models  described  in  Section  III  together  with  the  Forest  Fire 
reaction  rate  model  give  good  agreement  between  calculated 
velocity-loss  and  detonation  threshold  data  for  PBX-9404. 

For  Composition  B it  is  necessary  to  increase  the  Forest  Fire 
reaction  rates  by  a factor  of  1.5  in  order  to  obtain  agree- 
ment between  the  calculated  and  experimental  detonation  thres- 
holds. With  regard  to  the  factor  of  1.5,  since  the  Forest 
Fire  reaction  rates  were  obtained  with  the  HOM  equation  of 
state  for  partially  reacted  explosive  (Reference  4)  instead 
of  the  models  used  in  the  HELP  code  calculations,  a study  was 
undertaken  to  determine  the  sensitivity  of  the  reaction  rates 
to  the  particular  equation  of  state  models  used  in  their 
calculation.  The  results  of  this  investigation  are  given  in 
Section  4.2. 

2.2  SENSITIVITY  TO  UNREACTED  EXPLOSIVE  YIELD  STRENGTH 

Two  calculations  for  Composition  B were  performed  in 
order  to  determine  the  sensitivity  of  the  calculated  detona- 
tion threshold  to  the  value  of  the  yield  strength  used  for  the 
unreacted  explosive.  For  one  calculation  the  yield  strength 
was  increased  by  a factor  of  10  over  the  nominal  value  of 
0.512  kbar  (Section  3.2.1),  and  for  the  other  calculation 
the  yield  strength  was  reduced  to  zero.  There  was  no  dis- 
cernible effect  on  the  detonation  threshold  and  a negligible, 
if  any,  effect  on  the  calculated  velocity  loss.  It  should  be 
emphasized  that  this  insensitivity  to  the  value  of  the  yield 
stress  is  for  an  unconfined  explosive  geometry  at  a rather 
high  striking  velocity.  For  a cased  configuration  at  lower 
impact  velocities,  the  value  of  the  yield  strength  and  its 
effect  on  the  amount  of  plastic  deformation  would  probably 
be  significant. 
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III.  MATERIAL  MODELS 


3.1  INERT  MATERIAL  MODEL  AND  PARAMETERS  FOR  STEEL 


The  inert  material  model  employed  in  HELP  includes  an 
equation  of  state,  a deviatoric  constitutive  relation  for 
elastic  and  plastic  deformations,  and  a failure  criterion,  as 
discussed  in  the  following  paragraphs. 

For  steel,  the  Tillotson  (Reference  5)  equation  of 

state,  modified  to  give  a smooth  transition  between  compressed 

and  expanded  states  was  used.  For  compressed  states,  i.e., 

when  p/p«  > 1,  or  for  any  cold  states,  E < E , where  E is 
w s s 

the  specific  internal  energy  at  incipient  vaporization,  the 
equation  of  state  has  the  form 


P = P = 
c 


a + 


+ 1 


Ep  + Ari  + Bn 


(3.1) 


where  C = p/Pq  n = S - 1*  For  expanded  hot  states,  i.e., 
when  p/Pq  < 1 and  E > E^,  the  equation  of  state  has  the  form 


A smooth  transition  between  the  compressed  and  expanded 
states  is  insured  by  a transition  equation  for  the  inter- 
mediate region  defined  by  E < E < E'  and  p/Pf,  < 1,  where 

5 S U 

E'  is  the  specific  internal  energy  at  complete  vaporization. 
This  blended  portion  of  the  equation  of  state  has  the  form 
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P 


(e-e^)Pe  -h  (e;-e)p^ 


(3.3) 


The  values  of  the  parameters  as  used  for  steel  in  the  above 
equations  are 


A = 1280  kbar 
B = 650  kbar 
a =0.5 
b =1.5 

Eq  = 9.5  X 10^®  erg/g 

Eg  = 2.44  X 10^^  erg/g 

E!  = 10.2  X 10^°  erg/g 
s 

o =5.0 
6 = 5.0 


The  deviatoric  stress  increments,  are  calculated 

from  the  elastic  relation 


2wde ! . , 
13 


(3.4) 


where  v is  the  shear  modulus  and  de|j  are  the  increments  of 
deviatoric  strain.  When  such  an  increment  of  stress  causes 
the  von  Mises  yield  criterion. 


Sf  . < 2Y^ 
1]  - s 


to  be  violated,  each  stress  component  is  proportionately  reduced 
to  bring  the  stress  back  (normally)  to  the  yield  surface.  A 
variable  yield  strength 


R-3512 


= (vo  * V,.  * (l  - ^ ) 


(3.5) 


can  be  used  to  account  for  the  increase  in  strength  at  high 
pressures  and  the  decrease  in  strength  at  elevated  values  of 
the  specific  internal  energy,  where  E is  the  internal  energy 
at  the  melt  point.  The  values  of  the  parameters  as  used  for 
steel  are 


= 800  kbar 


= 6.0  kbar 


= 1.3  X 10  erg/g 


The  failure  criterion  employed  in  the  HELP  code  is 

based  upon  a minimum  allowable  compression  ratio,  n . . 

inxxi 

When  the  material  in  a computational  cell  is  subjected  to 
conditions  such  that  the  failure  criterion  is  violated,  i.e.. 


min 


the  material  is  considered  failed,  and  the  stresses  in  that 
cell  are  set  to  zero.  For  steel,  the  value  used  for  the 
minimum  compression  ratio  was 


ri  . = -0.02 

min 
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3.2  HIGH -EX PLOSIVE  EQUATION  OF  STATE 

Partially  reacted  explosive  is  treated  as  a mixture 
of  completely  reacted  and  completely  unreacted  components  in 
pressure  equilibrixim.  Separate  equations  of  state  of  the 
general  form 


P 


2 


f 2 ( p2 ' ^2 ^ 


(3.6) 


are  used  for  the  unreacted  and  reacted  components,  respectively, 
where  p2  and  E^^,  E2  are  the  densities  and  specific  internal 
energies.  The  equations  of  state  are  described  in  detail  in 
the  following  sections. 

The  pressure  is  calculated  as  follows:  at  the  end  of 
a computational  cycle  the  overall  density  p and  the  specific 
internal  energies  of  the  unreacted  and  reacted  components 
are  known.  The  general  expression  relating  the  overall  density 
to  the  individual  densities  of  the  unreacted  and  reacted  com- 
ponents is 


1 

"p 


(1-w) 


(3.7) 


where  w is  the  mass  fraction  of  unreacted  explosive  in  the 
computational  cell.  The  individual  densities  are  varied  in 
an  iterative  procedure,  subject  to  the  above  constraint,  until 
the  calculated  pressures  of  the  unreacted  and  reacted  components 
are  equal.  Thus,  the  determination  of  the  pressure  in  a cell 
containing  partially  reacted  high-explosive  is  formally  the 
same  as  that  for  a cell  containing  a mixture  of  any  other  two 
material  constituents,  and  the  computational  machinery  developed 
previously  for  use  in  the  HELP  code  was  immediately  applicable. 
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The  determination  of  the  mass  fraction  of  unreacted 
explosive  in  a computational  cell  is  a two-step  procedure. 

The  first  possible  variation  is  due  to  the  mass  flow  into  and 
out  of  a given  cell  from  neighboring  cells  with  different 
unreacted  mass  fractions.  This  increment  is  calculated  at  the 
same  time  that  the  cell  mass  is  updated,  a mass-weighted  aver- 
age of  the  appropriate  contributions  being  used.  The  second 
variation  is  explicit  change  as  calculated  from  the  assumed 
reaction  rate  of  the  explosive  decomposition.  A separate 
subroutine  was  written  to  perform  this  part  of  the  update. 

3.2.1  Unreacted  Composition  B 

For  compressed  states,  p > Pq,  a Mie-Gruneisen  equation 
of  state  was  used, 

= f(p)  + pGE,  (3.8) 


where 


f(p)  - (l  - f 5c)  P„(P) 


(3.9) 


in  terms  of  the  Hugoniot  pressure  P„(p),  and 

n 


r = £_  e = 1 - iL 
E - p ' ' - 1 Vj  • 


(3.10) 


The  Hugoniot  pressure  is  determined  from  the  straight  line 

U^,  U fit, 
s p 


= C,  + SUp  , 


(3.11) 
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or  equivalently, 

Ph  ' PqCq  e/(l  - Se)^  . (3.12) 

For  negative  pressures  the  Hugoniot  was  extrapolated  as  a 
straight  line.  Values  of  the  parameters  in  the  above  equa- 
tions as  obtained  from  Reference  1 are 

3 

Pq  = 1.715  g/cm 
Cq  = 2.31  X 10^  cm/sec 
S = 1.83 
G = 1.5  . 

Other  parameters  as  obtained  from  the  Ballistic  Research 
Laboratories  (Reference  6)  are  the  shear  modulus, 

V = 50.3  kbar 

the  yield  stress  in  simple  tension  of 
= 0.887  kbar 

from  which  the  von  Mises  yield  criterion  is 
Yg  « Y^//3  * 0.512  kbar, 

and  the  tensile  strength  of 
-S*  = 0.1042  kbar. 

The  volumetric  strain  at  failure  as  used  in  the  HELP  code  is 


'r' 

1] 
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t , 


; ! 


i I 
' ] 


i 


Li 


then 


^n.in  ' SVIP^CJ) 


- 0.00114 


For  expanded  hot  states  the  equation  of  state  used 


was 


Pj.  = OGE  4 - ij  exp  j-6  - l) 


exp  [-p(^-  l)  ] • 


(3.13) 


For  intermediate  states,  the  blend  as  described  in  Section 
3.1  was  used.  For  Composition  B,  the  energy  at  the  melt  point 
was  taken  as 


.8 


E = 8.29  X 10  erg/g 
m 


and  the  energy  at  incipient  vaporization  and  complete  vaporiza- 
tion were  estimated  as 


E = 9.19  X 10’  erg/g 
s 


ao 


E'  « 3.0  X 10  erg/g 
s 


The  standard  values  a “ 6 “ 5.0  were  used  in  Equation  (3.13) 


3.2.2  Unreacted  PBX-9404 

The  Mie-Gruneisen  equation  of  state,  as  described  in 
the  preceding  section,  was  also  used  for  compressed  states 
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of  PBX-9404.  Values  of  the  parameters  as  obtained  from 
Reference  7 are 

Pq  = 1.83  g/cm^ 

Cq  = 2.423  X 10^  cm/sec 

S = 1.883 
G = 0.675 

The  shear  modulus,  as  calculated  from  the  shear  wave  velocity 
of  1.62  X 10^  cm/sec  (Reference  8)  is 

2 

U = ^o^t  ” 48.0  kbar. 

The  form  of  the  expanded  phase  equation  of  state  was 
taken  to  be  the  same  as  that  for  Composition  B,  Equation  (3.13), 
with  a = 8 = 5.0.  The  melt  temperature  was  taken  as  280®C, 
which  together  with  specific  heat  data  from  Reference  9 gives 
for  the  energy  at  melt, 

E = 3.21  X 10®  erg/g. 
m 

The  energies  for  incipient  vaporization  and  complete  vaporiza- 
tion were  estimated  as 


E “ 1.16  X 10^®  erg/g 
8 
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The  yield  and  failure  parameters  were  taken  to  be  the 
same  as  those  for  Composition  B,  i.e.,  the  shear  yield  stress 
for  use  in  the  von  Mises  yield  criterion  is 

Yg  = 0.512  kbar, 

and  the  tensile  strength  as 
-S*  = 0.1042  kbar. 

The  resulting  volumetric  strain  at  failure  is  then 

"min  = SV(PoC^)  . - 0.00097 

3.2.3  Reacted  Composition  B 

The  Jones-Wilkins-Lee  (JWL)  equation  of  state  was  used 
for  both  Composition  B and  PBX-9404  for  the  completely  reacted 
state.  It  had  already  been  coded  for  use  in  the  HELP  code 
and  is  of  the  form 

^ ~ * (^  " R^)  + B ^1  - exp  (-R2V) 


where  V is  the  volume  of  the  detonation  products  relative  to 
that  of  the  undetonated  explosive,  and  E is  the  internal 
energy  per  cm^  of  undetonated  explosive. 

The  values  of  the  parameters  for  Composition  B as 
taken  from  Reference  10  are 

A « 5.242  Mbar 
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3.2.4  Reacted  PBX-9404  j 

i 

The  pareuneters  for  the  JWL  equation  of  state  as  obtained  | 

from  Reference  10  are 

A =8.545  Mbar  | 

B = 0.2049  Mbar  I 

j 

R.  = 4.60  1 

^ \ 

R^  = 1.35  I 

to  = 0.25  j 

Eq  = 0.102  Mbar-cc/cc 


f 
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IV.  REACTION-RATE  MODEL 


FOREST  FIRE  REACTION-RATE  MODEL 


A brief  description  of  the  model  is  given  here;  a de- 
tailed discussion  is  given  by  Mader  and  Forest  (Reference  1) . 

The  basic  assumption  is  that  the  build-up  to  complete 
detonation  takes  place  along  a unique  curve  in  distance,  time 
and  state  space,  regardless  of  the  initial  driving  pressure. 
The  shock  pressure  as  a function  of  distance  run  is  then 
given  by  the  so-called  "Pop  plot",  which  is  a fit  to  the  ex- 
perimental data  for  initial  shock  pressure  as  a function  of 
run  to  detonation.  The  Pop  plot  for  Composition  B is  shown 
in  Figure  4 along  with  the  assumed  growth  of  the  shock  pres- 
sure as  complete  detonation  is  approached.  A second  assump- 
tion, as  indicated  in  Figure  4,  is  that  the  pressure  gradient 
immediately  behind  the  shock  front  is  zero. 

The  additional  information  required  for  the  determina- 
tion of  the  reaction  rates  are  the  Hugoniot  curve  for  shocks 
to  the  intermediate,  partially  reacted  states,  and  an  equation 
of  state  of  the  form 


P — f (p,E,w), 


(4.1) 


where  p is  the  density,  E the  specific  internal  energy,  w the 
mass  fraction  of  unreacted  explosive.  The  calculational  pro- 
cedure is  as  follows: 

For  a given  shock  pressure  P^j,  the  experimental  curve 

of  U versus  U for  the  partially  reacted  Hugoniot, 
s p 


U « C + S'U  , 
s p 


(4.2) 


and  the  Hugoniot  relations  are  used  to  determine  the  state 
at  the  shock  front, 


“ ^oVp  - PQ  * S*  Up) Up 


(4.3) 
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V„  = V (U  - U^)/U 

n Os  p S 


E = i u2 

2 • 


(4.4) 


(4.5) 


These  values  are  used  with  the  equation  of  state  to  determine 
the  mass  fraction  of  unreacted  explosive  at  the  shock 
front , 


(4.6) 


The  reaction  rate  itself  is  then  obtained  from  the  Lagrangian 
time  derivative  of  the  above  expression, 


(i ) = (If)  * (lf)„  * ^3  (If)  • 


(4.7) 


where  the  subscript  m denotes  the  time  derivative  at  a fixed 
mass  element,  and  f.  _ _ are  the  partial  derivatives  of  the 

± f /.  f O 

equation  of  state  function  with  respect  to  its  arguments. 

At  this  point,  the  assumption  of  a zero  pressure 
gradient  at  the  shock  front  is  used  to  evaluate  the  partial 
derivatives.  The  general  expression 


©.  ■ (SI.  • (S).  (g), 


reduces  to 


/9P\  ^ ^ ^ 

[Btj^  dt  " d1 


(4.8) 
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where  the  derivatives  of  the  Hugoniot  pressure  with  respect 
to  R (run  to  detonation)  is  obtained  from  the  Pop  plot.  The 
Lagrangian  equation  for  conservation  of  momentum 


(’t)m  ■ ' (*”)t  ‘ “ 

allows  one  to  write  for  the  rate  of  change  of  particle 
velocity  at  the  shock  front 

= /9u\  . /3u\  ^ „ p /au\ 

dt  = (at)^  (am)^  dt  - pQ^s  (amj^  ' 


(4.9) 


and  the  Lagrangian  equation  for  conservation  of  mass  yields 


at 


(4.10 


for  use  in  Equation  (4.7),  where  the  time  derivative  on  the 
right  in  the  above  expression  is  evaluated  from  the  pressure 
time  derivative  by 


^ = Vo3^/(s'Uptu^). 


(4.11 


The  energy  derivative  in  Equation  (4.7)  is  obtained  from  the 
Lagrangian  equation  for  conservation  of  energy, 
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4.2  SENSITIVITY  OF  REACTION  RATE  TO  EQUATION  OF  STATE 

Since  the  Forest  Fire  reaction  rates  were  obtained  us- 
ing the  HOM  equation  of  state  (Reference  4)  in  Equation  (4.1) 
to  evaluate  the  pressure  for  a partially  reacted  state,  it 
seemed  worthwhile  to  recompute  them  using  the  HELP  models  in 
order  to  check  the  sensitivity  to  variations  in  the  explo- 
sive equation  of  state.  Such  a calculation  was  done  for 
Composition  B. 

As  a first  step,  the  Fortran  lisings  given  in  Reference 
1 were  punched,  and  the  program  was  run  using  the  original  HOM 
equation  of  state  in  order  to  insure  that  there  were  no  coding 
errors.  The  reaction  rates  for  Composition  B were  then  re- 
calculated using  the  equation  of  state  models  described  in 
Section  III  and  the  appropriate  subroutines  from  the  HELP 
code.  Comparison  of  Hugoniots  for  several  unreacted  mass 
fractions  is  given  in  Figure  5,  and  a comparison  of  the 
reaction  rates  as  calculated  from  the  two  equations  of  state 
is  given  in  Figure  6.  As  the  figure  shows,  the  differences 
of  the  equations  of  state  are  enough  to  very  nearly  account 
for  the  factor  of  1.5  necessary  to  bring  the  experimental 
and  calculated  detonation  thresholds  for  Composition  B into 
agreement. 
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Hugoniots  for  Composition  B 


HOM  Equation  of  State 


HELP  Equation  of  State 
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Figure  5. 


Volume  (cc/g) 


A comparison  of  several  Hugoniots  for  Composition 
B as  given  by  the  HOM  equation  of  state  and  the 
models  used  in  the  HELP  code  calculation,  where  w 
is  the  mass  fraction  of  unreacted  explosive. 


Reaction  Rate  Msec 
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V.  CONCLUSIONS 

For  the  explosives  PBX-9404  and  Composition  B the  use 
of  the  equation-of-state  models  presented  here  and  the  Forest 
Fire  reaction  rates  (modified  for  Composition  B to  account 
for  the  sensitivity  to  the  equation  of  state)  results  in  good 
agreement  between  calculated  and  observed  velocity  losses  and 
detonation  thresholds  for  the  penetration  of  unconfined  explo- 
sive slabs  by  steel  spheres.  The  reaction-rate  model  (Reference 
1) , which  was  published  after  the  start  of  the  present  contract, 
has  also  been  used  with  success  at  the  7<os  Alamos  Scientific 
Laboratory  to  describe  the  detonation  process  in  heterogeneous 
explosives . 

However,  as  pointed  out  in  Reference  1,  the  use  of  this 
reaction-rate  model  is  not  realistic  at  very  low  pressures. 

For  example,  in  the  calculations  reported  here,  the  reaction 
rates  are  assumed  to  be  a function  of  pressure  only,  and  they 
are  set  to  zero  for  pressures  less  than  15  kbar  in  PBX-9404 
and  20  kbar  in  Composition  B.  Consequently,  the  model  as 
described  here  would  not  be  useful  for  confined  geometries 
at  low  pressures,  in  which  the  explosive  has  a chance  to 
"cook" , or  for  configuration  in  which  the  initial  heating  is 
due  principally  to  plastic  deformation  associated  with  macro- 
acopic  shear . 
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